Effect of a magnetic field on the quasiparticle recombination in superconductors 
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Quasiparticle recombination in a superconductor with an j-wave gap is typically dominated by a phonon 
bottleneck effect. We have studied how a magnetic field changes this recombination process in metallic 
thin-film superconductors, finding that the quasiparticle recombination process is significantly slowed as 
the field increases. While we observe this for all field orientations, we focus here on the results for a 
field applied parallel to the thin film surface, minimizing the influence of vortices. The magnetic field 
disrupts the time-reversal symmetry of the pairs, giving them a finite lifetime and decreasing the energy 
gap. The field could also polarize the quasiparticle spins, producing different populations of spin-up and 
spin-down quasiparticles. Both processes favor slower recombination; in our materials we conclude that 
strong spin-orbit scattering reduces the spin polarization, leaving the field-induced gap reduction as the 
dominant effect and accounting quantitatively for the observed recombination rate reduction. 
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PACS numbers: 74.40.Gh, 74.70.Ad, 78.47.D-, 74.25.Ha 

An excitation from the superconducting condensate re- 
quires finite energy (the energy gap 2A) and produces two 
quasiparticles. A quasiparticle excited to very high energy 
(compared to A) quickly decays via a number of fast scat- 
tering processes to near the gap edge, where it recombines 
with a partner to form a Cooper pair. Thepair's binding 
energy is emitted mainly as 2 A phonons [Ll] — [3D - The re- 
combination process is delayed by a phonon bottleneck: 
each recombination-generated phonon can break another 
Cooper pair, causing energy to be trapped in a coupled 
system of 2A phonons and excess gap-edge quasiparti- 
cles lll|5llal- Quasiparticle recombination has been widely 
studied in both metallic superconductors, to investigate the 
non-equilibrium processes in the many-body BCS system 
ffvl— TToh . and high-temperature superconductors, to gain new 
insight into the pairing mechanism |[llT - [l4ll . Theories of 
the recombination process considered the reaction kinet- 
ics and interactions of quasiparticles and phonons 
while experiments obtained the dependence of the quasi- 
particle lifetime on temperature, film thickness, and excita- 
tion strength IU7U18II . A magnetic field is known to couple 
to the electron orbital motion and to align the spin; both ef- 
fects weaken superconductivity 111 911 . The consequence of 
magnetic field on the quasiparticle recombination 12011 has 
not been examined in detail by optical pump-probe meth- 
ods. 

We use a novel time-resolved laser-pump synchrotron- 
probe spectroscopic technique to study the quasiparticle re- 
combination dynamics in superconducting thin films, un- 
der applied magnetic field. Samples studied include a 
10 nm thick Nbo.5Tio.5N film on a crystal quartz substrate 
and a 70 nm thick NbN film on a MgO substrate. These 
substrates are essentially transparent in the far-infrared 
spectral range. The films were grown by reactive mag- 
netron sputtering, using NbTi cathode in Ar/N 2 gas for 
Nbo.5Tio.5N and Nb cathode in N 2 gas for NbN. The two 



films have critical temperatures of 10.2 K and 12.8 K, and 
a zero-temperature, zero-field, energy gap 2A of 2.7 meV 
and 4.5 meV, respectively. Four-probe resistivity measure- 
ments with magnetic field parallel to the films determined 
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FIG. 1. Experimental set-up. Electrons circulate in bunches in 
the synchrotron storage ring, generating pulses of far-infrared ra- 
diation with a repetition frequency of 52.9 MHz. The Ti:sapphire 
laser produces pulses with a repetition frequency of 105.8 MHz 
and a pulse picker selects every other pulse to match the syn- 
chrotron pulse pattern. The selected laser pulses are delivered 
over a fiber optic cable to the sample and the synchrotron pulse 
probes the photoinduced transmission at a fixed time delay af- 
terward. To synchronize the synchrotron and laser pulses, the 
52.9 MHz bunch timing signal from a pair of electrodes inside 
the synchrotron ring chamber is used by the Synchro-Lock laser 
control system as a reference for the laser pulse emission. The 
pulse generator introduces an adjustable delay between the laser 
and synchrotron pulses. The transmitted far-infrared light is de- 
tected by a bolometer detector and recorded on a computer. 
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FIG. 2. Photoinduced transmission S(t) vs. time t for 
Nbo.5Tio.5N (a and b) and for NbN (c and d), all measured in 
parallel fields at T < 2 K. Low-fluence and high-fluence data 
are compared. Note the semilog scale; simple exponential decay 
produces a straight line. 
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their upper critical field to be greater than 20 T at T < 3 K. 

The samples were mounted in a 4 He Oxford cryostat 
equipped with a 10 T superconducting magnet, and probed 
by far-infrared radiation produced in a bending magnet at 
beamline U4IR of the National Synchrotron Light Source, 
Brookhaven National Laboratory. The experiment, illus- 
trated in FIG. [H exploits the fact that the synchrotron ra- 
diation is emitted in ~300 ps long pulses (governed by the 
electron bunch structure in the storage ring). We applied 
mode-locked near-infrared Ti:sapphire laser pulses (~2 ps 
in duration and ~1.5 eV in photon energy) as the source 
for photoexcitation. The synchrotron probe beam measures 
the photoinduced optical properties due to the excess quasi- 
particles as a function of time delay relative to the arrival 
of the pump beam. The synchrotron pulse has a Gaussian 
profile with a FWHM of ~300 ps, determining the time 
resolution of the experiment. At selected delay times t, we 
measure the spectrally integrated photoinduced transmis- 
sion S(t) = —AT(t) over the spectral range spanning the 
superconductors energy gap (~3 meV). The spectral shape 
is determined primarily by the optical components carrying 
the beam, and the detector. 

If the laser were turned on and off to measure the pho- 
toinduced response, there would be a temperature modu- 
lation as well as the photoexcited quasiparticle modula- 
tion. To reduce these thermal effects we dither the laser 
pulse back and forth by a few tens of picoseconds at each 
delay setting, keeping the incident laser power constant. 
The dither is achieved by phase modulating the laser pulse 
using the internal oscillator of a lock-in amplifier. The 
directly obtained quantity is therefore a differential sig- 
nal, dS/dt. This signal was detected using a B-doped Si 
bolometer in combination with the lock-in amplifier. Nu- 
merical integration yields the photoinduced transmission 
S(t), which directly follows the excess quasiparticle den- 
sity 

To study the effect of magnetic fields and excess carrier 



FIG. 3. Effective instantaneous recombination rate vs. pho- 
toinduced transmission, (a) For Nbo.5Tio.5N, data at each field 
include fluences ranging from 0.4 to 10.7 nj/cm 2 . (b) For 
NbN, data at each field include fluences ranging from 2.4 to 
18.1 nj/cm 2 except for 8 T and 10 T, where data were collected 
at 18.1 nj/cm 2 . A 4-point moving average was performed on the 
data to reduce noise. The lines are linear fits to the data. 

density on the recombination dynamics, we measured the 
magnetic -field and laser-fluence dependent photoinduced 
transmission for Nbo.5Tio.5N and NbN thin films. The sam- 
ples were fully immersed in superfluid 4 He (T < 2 K) 
to minimize heating. At this low temperature, the thermal 
quasiparticle population is small (compared to the number 
of broken pairs at high fluence) but not zero. The field 
was applied parallel to the film surface to avoid the com- 
plexity of vortex effects. (See Ref. yfl). Typical results 
are shown in FIG.|2] where the photoinduced signal S(t) 
(excess quasiparticle density) is plotted against delay time. 
At both low (FIG.|2k and FIG. [2]:) and high laser fluences 
(FIG. I2J3 and FIG. |2ji), a longer time is required for re- 
combination as the magnetic field is increased. The pulse 
width of the synchrotron probe beam gives rise to the initial 
upturn in the data, which is skipped in the following data 
analysis. 

We have discovered a revealing perspective to dis- 
play our results, shown in FIG. [3j We define an ef- 
fective instantaneous recombination rate l/r e g(t) = 
— [dS(t)/dt]/S(t) and plot l/r e s(t) vs S(t) at various 
fields and fluences. Here short times are at the right (large 
S(t)) and long times at the left. In this presentation, data at 
the same field but for different pump fluences scale to the 
same straight line. As will be shown below, this behavior 
is expected for bimolecular recombination where the life- 
time for a given particle is proportional to the availability 
of other particles with which to combine. 

The scaling can be understood as follows. The phonon 
bottleneck was first discussed by Rothwarf and Taylor 0]] 
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FIG. 4. Panels a and b show the excitation gap, (squares) and the pair-correlation gap A (triangles) for Nbo.5Tio.5N and NbN, 
obtained from the optical conductivity (left scale). The solid lines are theoretical calculations of A and Qq. The square root of the 
condensate density \ZN SC (proportional to the order parameter) is shown as circles (right scale). Panels c and d show the slope extracted 
from FIG.[3]vs. Qq from a and b. The error bars in both plots are calculated deviations of the slope from the linear fit in FIG. [3] The 
lines are linear fits to the circles. 



using two rate equations, one for the quasiparticles and 
the other for the 2A phonons. The quasiparticles, which 
directly correspond to our signal S(t), follow a simple 
model that captures the feature of bimolecular recombi- 
nation, meanwhile taking into account the phonon bottle- 
neck. The decay rate of the total quasiparticle density N(t) 
toward the equilibrium density is proportional to N 2 , be- 
cause recombination requires the presence of two quasipar- 
ticles. Motivated by the Rothwarf -Taylor yj] equations, we 
write 



d -^ = -2R{N 2 -Nl). 



(1) 



A thermal term iV t 2 h is subtracted from N 2 , because at equi- 
librium N = N til and the quasiparticle density must re- 
main constant. The phonon bottleneck is introduced into 
the model through the recombination rate coefficient R. 
(See Section 1 of the Supplemental Material.) A factor 
of 2 is included because each recombination event de- 
pletes two quasiparticles. Now, N(t) = N th + N ox (t), 
with N th the thermal density and N ex (t) the photoin- 
duced excess density. At a given temperature and magnetic 
field, N th is time -independent, making Eq. {D become 
-{dN D Jdt)/N cx = 2R(N cx {t) + 2iV th ). We identify 
— (dN cx / dt) /iV ex: as the effective instantaneous relaxation 
rate 1/T e g(t) defined earlier, because the photoinduced 
transmission S(t) is proportional to the excess quasipar- 
ticle density [2], S = CN CX , where C is just a constant to 
convert from signal to quasiparticle density. Hence, 

» «9M-» (SW + SSWl (2, 



S(t) dt C 
with 5 t h = CN th . Eq. (O is consistent with the linear 



behavior demonstrated in FIG. [3j The field dependence 
requires the prefactor R to decrease with field. 

To interpret the field dependence shown in FIG. [3] it 
is a prerequisite to understand how the field changes the 
electronic states of the superconductor. If spin-orbit scat- 
tering is small, the magnetic field could make the major- 
ity of quasiparticles have one spin direction. (This is the 
same polarization that gives Pauli paramagnetism to met- 
als.) Spin polarization will slow the recombination because 
only quasiparticles with opposite spins can recombine. A 
recombination model including this spin polarization ef- 
fect is discussed in Section 2 of the Supplemental Mate- 
rial. In this case, the recombination equation remains in 
the same form as Eq. ©, but with the coefficient 2R/C 
replaced by (8R/C)(N^ N± /N 2 ), where JV* and are 
respectively the densities of spin-up and spin-down quasi- 
particles. The quasiparticle spin-polarization factor N^N^ 
would depend on the magnetic field in the limit of weak 
spin-orbit coupling, just as in the Pauli paramagnetism of 
metals. According to the BCS theory, electrons form spin- 
singlet pairs condensed in the ground state; the spin suscep- 
tibility vanishes as the temperature approaches 0. The stud- 
ies of superconductor spin susceptibility were done on thin 
films with thickness so small that the effect of a magnetic 
field on the electron orbit could be neglected. Paramag- 
netic splitting of the quasiparticle density of states was ob- 
served in 5 nm aluminum films in a parallel magnetic field 
Ir23ll . In a study of magnetic field effects on far-infrared ab- 
sorption of thin superconducting aluminum films, van Ben- 
tum and Wyder [24] concluded that paramagnetic splitting 
was important in their thinnest films, but did not allow for 
quasiparticle spin polarization. If a high degree of spin po- 
larization existed, the recombination rate would be slowed 
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much more than observed. However, spin-orbit scatter- 
ing must be considered. Tedrow and Meservey observed 
the spin-state mixing in thin aluminum films due to spin- 
orbit scattering [25]. They defined a spin-orbit scattering 
parameter b = H/3At so to describe the degree of spin- 
orbit scattering, where r so is the spin-orbit scattering time. 
They calculated that, as b is increased to 0.5, the spin-up 
and spin-down quasiparticle density of states completely 
mix, leaving no clear signature of the two-peak feature in 
the density of states due to Zeeman splitting. Considering 
the short spin-orbit scattering time measured II26I1 in NbTi, 
r so = 3.0 x 10~ 14 s, and using the A of Nb 5 Ti 5 N and 
NbN, we estimate that b = 4.2 and 3.3 for Nb . 5 Ti . 5 N 
and NbN, respectively. We believe that spin is not a good 
quantum number in our samples, requiring us to look be- 
yond spin polarization to understand the recombination. 

In a study yfl of the optical conductivity of Nbo.5Tio.5N, 
we found that a parallel magnetic field breaks the time- 
reversal symmetry of the Cooper pairs and decreases the 
superconducting energy gap. The physics is similar to 
magnetic-impurity-induced pair-breaking effects, as origi- 
nally formulated by Abrikosov and Gorkov Il27ll . In a mag- 
netic field, one must distinguish between the spectroscopic 
energy gap 2f2 G and the pair-correlation gap A. These 
gaps 14J] are plotted in FIG. 0k as squares and triangles re- 
spectively. The real part of the optical conductivity, corre- 
sponding to the electromagnetic absorption, shows a clear 
suppression of the energy gap 2Q G with field (squares in 
FIG. The imaginary conductivity is a measure of the 
superconducting condensate density N sc , which goes as 
A 2 . The field dependences of A and of y/N sc (shown as 
circles) agree well, providing clear evidence for a weaken- 
ing of superconductivity by the magnetic field. The quan- 
tities S7 G , A, and N sc for NbN, obtained using the same 
technique (in Section 3 of the Supplemental Material), are 
plotted in FIG. @J}. The NbN field dependence is quali- 
tatively different from that of Nbo.5Tio.5N because in this 
thicker film the applied field induces a spatial variation 
in the order parameter, making the weakening of super- 
conductivity be proportional to the field, rather than being 
quadratic in field as in the much thinner Nbo.5Tio.5N lE9ll . 
The energy gaps will be used in the following analysis. 

The field dependence, shown in FIG.[3j is dominated by 
the recombination rate coefficient R. On the one hand, by 
explicitly solving Eq. (0 one can identify a low-fluence 
recombination rate l/r e g = ARN th . (See Section 4 of 
the Supplemental Material.) The field dependence of the 
thermal quasiparticle density N th results from the field- 
dependent energy gap and the quasiparticle density of 
states [30]. On the other hand, the effective lifetime of the 
excess quasiparticles is modified from the intrinsic value 
Tr, and is tied to the rates at which the phonons, produced 
in recombination events, re -break pairs (1/tb) or escape 
from the film (l/r 7 ) @]. The quasi-equilibrium values of 
Tr and t b were derived by Kaplan et al. [6]. Magnetic - 
field-induced pair breaking decreases the spectroscopic en- 



ergy gap (FIG. 0]l and modifies the quasiparticle density 
of states, resulting in a decrease in tr and an increase in 
Tb- The field independent phonon escape time is deter- 
mined by the film thickness and the acoustic mismatch be- 
tween the film and the environment lf32h . The recombina- 
tion rate coefficient R (and, hence, the slope of 1 /r e g in 
FIG. [3| is therefore field-dependent through 7V th , Tr and 
Tb- (See FIG. S7 in the Supplemental Material.) The equa- 
tion is involved but, when we compute the slope vs Q G for 
Nbo.5Tio.5N and NbN, shown in FIG.gh and FIG. Si, we 
obtain a basically linear relation. This calculation implies 
a connection between the field-dependent quasiparticle re- 
combination and the field-induced pair breaking. The lin- 
ear relation can be explained by considering only the field- 
induced gap reduction. (See FIG. S7 in the Supplemental 
Material.) The finite y-intercepts in FIG. Hfc and FIG. EJl 
are intriguing, bringing out the question of how the pho- 
toexcited quasiparticles relax in a gapless superconductor, 
motivating challenging experiments to probe the gapless 
regime. 

In conclusion, our time-resolved pump-probe measure- 
ments on metallic s-wave superconductors reveal a slow- 
ing of quasiparticle recombination in an external magnetic 
field. The field was aligned parallel to the thin-film sam- 
ple surface, to minimize effects due to vortices. There are 
two possible causes of the observed slowing: field-induced 
spin imbalance and field-induced gap reduction. The spin 
imbalance is unlikely to be important in Nb 5 Ti 5 N and 
NbN due to strong spin-orbit scattering. This scenario can 
be tested by investigating materials with small spin-orbit 
scattering. The field-induced gap reduction alone can ex- 
plain quantitatively the slowing of recombination, and we 
conclude it to be the dominant effect observed in our ex- 
periment. 
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Supplemental Material 

1. Phonon bottleneck effect 

In this section we show that the phonon bottleneck effect 
can be included in the effective recombination rate R as a 
proportionality coefficient. Consider the Rothwarf-Taylor 
equations OJ] for the coupled populations of quasiparticles 
and phonons in the absence of external quasiparticle injec- 
tion, 

— = -Tk 1 N 2 + t b 1 N lu , (SI) 

^ = l^Ni-K^N^-T^-N^), (S2) 

where N(t) and N u (t) are respectively the densities of 
quasiparticles and high-energy (defined as hw > 2A) 
phonons, r^ 1 is the intrinsic quasiparticle recombination 
rate, r^ 1 is the phonon pair-breaking rate, and r^ 1 is the 
phonon escape rate. (The phonons may enter the sub- 
strate, enter the helium bath, or decay anharmonically 
to energies hw < 2A.) In thermal equilibrium, the 
quasiparticles N(t) and high-energy phonons N u (t) reach 
time-independent equilibrium values, linked by iV th = 

^ N^^Tg 1 /t^, where N th and iV^ th are for the quasi- 
particles and high-energy phonons, respectively. The cou- 
pled non-linear equations dSlI) and (IS2I) can be solved nu- 
merically, e.g., using the Runge-Kutta integration method. 
To illustrate the solutions, we set r^ 1 = 1 ns -1 and Tg 1 = 
10 ns -1 , typical values for these quantities at low temper- 
atures [2]. Without loss of generality we set N u<th = 1, 
which determines N th = \/T0. The initial value of N is 
determined by the pump laser fluence. We consider three 
cases, AT(0) = 2N th for low fluence, N(0) = 5N th for 
intermediate fluence, and -/V(O) = 10iV th for high fluence. 
iV w (0) is set to iV^th. For each case, we consider a range 
of phonon escape rates (t^ 1 ) from 1 ns -1 to 10 ns -1 , span- 
ning the strong to the weak phonon-bottleneck regime. 

The numerical solutions of the Rothwarf-Taylor equa- 
tions are shown in the first row of FIG. [Si] After a short 
period of approximately 0.05 ns, the phonon bottleneck 
effect becomes clear. As expected, when the phonon es- 
cape rate increases, the phonon bottleneck effect becomes 
weaker and recombination becomes faster. In the second 
row of FIG. [Si] the rate of change of the excess quasipar- 
ticle density N cx (t) = N(t) - N th is plotted vs. iV cx . 
The early stage (t < 0.05 ns) shown in the first row of 
FIG. [Si] has been skipped, because the quasiparticle and 
phonon populations are not yet equilibrated and it does not 
give information about the phonon bottleneck. Moreover, 
this stage is not temporally resolved in our experiment. 

The quasi-linear relation shown in the second row of 
FIG. ED suggests we can define an effective recombina- 
tion rate R as — (dN ex /dt) /N cx . This rate, defined as the 



slope of —dN cx /dt vs. A cx , is plotted vs. r 7 1 in FIG.IS2I 
As A(0) increases from 2N th to lOA^, the relation be- 
tween R and r" 1 becomes almost linear. We estimated 
that, for our experimental conditions, A(0) in our sam- 
ples was three orders of magnitude greater than A th at the 
lowest laser fluence. Therefore it is safe to conclude that 
the effective recombination rate R scales linearly with the 
phonon escape rate r" 1 . 

2. Pauli paramagnetism 

Our recombination model can be extended to include the 
magnetic-field dependence due to quasiparticle spin po- 
larization. In a quasiparticle recombination event, both a 
spin-up and a spin-down quasiparticle are needed to form 
a Cooper pair. We use an equation similar to the band- 
to-band recombination equation in a semiconductor to de- 
scribe the quasiparticle recombination, 

dN dN^ dN± . + , t . , 

*- - -* + -* = ~ 8R{N A ' ~ «»• (S3) 

where j and J. denote the spin-up and spin-down pop- 
ulations, respectively. This reduces to Eq. (1) in the 
main text when the spin-up and spin-down populations are 
equal. At a given condition, Eq. (IS3I) can be rewritten 
as -(dN cx /dt)/N cx = 8R(N^N+)(N ex + 2N th )/N 2 , 
where N cx = N — N th . Relating the measured photoin- 
duced transmission S to the excess quasiparticle density, 
S = CN CX , where C is a constant, the recombination 
equation can be further reduced to 

"If = ^(iVW)(S + 2S, h ), ,S4) 

where S th = CN th . 

The majority spin fraction N^/N can be calculated from 
the paramagnetic model, assuming Fermi-Dirac distribu- 
tion of quasiparticles f(E) and the quasiparticle density of 
states D(E) from the BCS theory, 

N l = 2 / f(E)D(E + nH)dE, (S5) 
Jo 

poo 

A T = 2 / f(E)D(E - fjLH)dE. (S6) 
Jo 

If we assume that the field dependence we see in FIG. 3 in 
the main text is only through the product N^N^, we can 
extract the majority spin fraction at different fields. The 
results are compared with the calculation in FIG.IS3I Pure 
Pauli paramagnetism predicts a stronger magnetic field de- 
pendence than observed in our data. 

The theory, however, must consider the strong spin-orbit 
scattering in NbN and Nbo.5Tio.5N. In the main text, we es- 
timate the strong spin-state mixing in the Nbo.5Tio.5N and 
NbN samples. Based on that argument and the analysis 
shown in this section, we expect that the spin polarization 
factor N^N^ is weakly dependent on the field for our sam- 
ples. 
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FIG. SI. Numerical solutions of the Rothwarf-Taylor equations. The first row shows the solutions for the quasiparticle density N and 
the high-energy phonon density iV w at different fluences and for a range of phonon escape rate. The second row shows the rate of 
change of the excess quasiparticle density — (dNcx/dtf/Ncx as a function of A^x- 
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FIG. S2. The effective recombination rate as a function of the FIG. S3. Majority spin fraction. The theory is based on Eqs. dS5b 
phonon escape rate at low, intermediate and high fiuence. and l lS6l >. The data points are extracted from the slopes in FIG. 3 

in the main text, all scaled so that the majority spin fraction is 1/2 

at T. 

3. Analysis of optical conductivity for NbN 



We studied the magnetic-field-induced effects in the 
Nbo.5Tio.5N and NbN thin films using Fourier transform 
far-infrared spectroscopy. The experimental technique and 
the analysis for Nbo.5Tio.5N can be found in Ref. [Q]. 
Here we analyze the superconducting-state to normal-state 
transmission ratio T s /7~ n and reflection ratio !Z s /lZ n for 



NbN shown in FIG. [S4l measured at 2 K with the mag- 
netic field parallel to the film. The data were taken with 
4 cm -1 (0.5 meV) resolution, so that the fringes due to the 
multiple internal reflections in the substrate were not re- 
solved. The angle of incidence for both transmission and 
reflection was 30°. The NbN thin film has a normal-state 
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FIG. S4. The transmission and reflection of NbN in parallel fields 
and 2 K, normalized to the corresponding normal-state values. 
The angle of incidence 0, for both transmission and reflection 
was 30°. 



conductivity a n = 2.0 x 10 3 Ohm - cm -1 , determined 
from its normal-state transmittance and thickness. The 
MgO substrate has a refractive index n ~ 3.0 and neg- 
ligible absorption in the far-infrared. The zero-field gap 
A = 17.9 cm -1 is obtained from fitting the zero-field 
optical conductivity with the Mattis-Bardeen theory. 

From the transmission and reflection ratios, we extracted 
the real (o{) and imaginary (ct 2 ) parts of the optical con- 
ductivity using the method discussed in Ref. yfl. The re- 
sults are shown in FIG. [S3] We found that the field depen- 
dence can be explained well by the pair-breaking theory. 
(See details of the theory in Ref. fjjj.) 

The pair-breaking parameter T is the only fitting param- 
eter, describing the strength of pair breaking due to the 
magnetic field. Its value at different fields is shown in 
FIG.ISol From T we calculated the pair-correlation gap A 




Frequency [cm 1 ] 

FIG. S5. The real (circles) and imaginary (triangles) parts of the 
optical conductivity of NbN at various parallel fields and T — 
2 K, normalized to the normal-state conductivity a n . The solid 
lines are fits to ai/a n using the pair-breaking theory Q]. The 
dashed lines show the corresponding 02/ 'cr n as determined by a 
Kramers-Kronig transform of the fit to the real part. 
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FIG. S6. Field dependence of the pair-breaking parameter Y de- 
termined from the optical conductivity of NbN, normalized to the 
zero-field gap. The dashed line is a linear fit. 



and the effective spectroscopic gap f2 G using ln(A/A ) = 
-tiT/4A and Q G = A[l - (T/A) 2/3 ] 3/2 . These quan- 
tities are shown in FIG. 4b in the main text. The super- 
conducting condensate density N sc is estimated from the 
below-gap part of cr 2 at T <C T c , which has the form, 



02 (w) 



iVsce 2 



moo 



(S7) 



where e and m are the electron charge and mass, respec- 
tively. 



4. Exact solution of the recombination equation 

The recombination equation, Eq. (2) in the main text, 
links the measured S(t) to the model, 

where S t h = CN th . This equation has the following exact 
solution: 

2e-*/ r 



S{t) = 5(0) 



2 + »l-e- 



(S8) 



t/r) 



Here r is the effective lifetime: r = l/4RN th . In the 
low-fluence regime, especially when 5(0) -C 5th, the so- 
lution is close to an exponential decay. As the fluence in- 
creases, the deviation from a simple exponential decay be- 
comes significant. 



5. Effective recombination rate 

On the one hand, the quasiparticles, interacting with 
phonons in the system, decay with an effective lifetime 
T G ff = t 7 + (1/2)t h (1 + Tj/tb), where r 7 , r R , and r B 
are the same quantities as defined above and in the main 
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text J3]. At low temperatures, t r 3> r 7 and r 7 3> r B 
|0]. The effective lifetime can be approximated as r c ff = 
t r t^/2t b . On the other hand, by solving the recombina- 
tion equation proposed in the main text, one can identify an 
effective quasiparticle lifetime r Gf j = l/4RN th , as shown 
in the previous section. As a result, 



R 



th 



T~ X T~ X 

2r B 1 iV th ' 



(S9) 



The phonon escape rate r 7 1 is expected to be independent 
of the gap, as discussed in the main text. A theory for Tr 1 
and Tg 1 has been given by Kaplan et al. |@], 



n 2 (n 



1 + 



(A; B T C ) 3 [1 - f(uj)\ J u+A — u) 2 - A 2 
A 2 



0,ph 



[n(fi) + l]/(fi-w)dn, (S10) 
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ttA(O) J A 
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Vw 2 - A 2 
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where Tq" 1 and r ~p h are respectively the characteristic life- 
times of the quasiparticles and phonons, determined by 
the electron-phonon coupling and phonon density of states. 
The quantities / and n are the Fermi-Dirac and Bose- 
Einstein distribution functions, respectively, r^ 1 should 
be evaluated at the quasiparticle gap energy u = A and 
Tg 1 should be evaluated at the phonon energy = 2 A for 
the estimation of their near-equilibrium values. For t r 1 

-1 /-oo Q3 p -(n-A)/k B T 

V(A) ' - - 



(k B T c 



dU- , 

,a ^(n- a) 2 



'0 e 



(a + 2A) 



- A 2 

5/2 



^-x/k B T 



(S12) 



in which we have replaced the Bose factor n(Q) and the 
Fermi factor /(A) by their low-temperature values. Be- 
cause only small x in the integrand contributes significantly 
to the integral, an expansion of the numerator yields 
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FIG. S7. The dependence of tr, tb, N^, and R, on the spec- 
troscopic gap, calculated directly from Eqs. dS 10b . dS lib . dS 15b . 
and dS9b without approximations. 



The phonon pair-breaking rate t b 1 has a simple form for 
its near-equilibrium state, given in Ref. I6J] as 



■ X (2A) 



A 



[l-2/(A)]«r, 



'0, P h A 

The quasiparticle density is 

N th « A^(0)V27rAfe B re- A/fcBT 
Substituting Eqs. (fST3T) — (fST3T > into Eq. 
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2^2A r-Vo 1 1 
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5k B T 
8A 



A. 
(S16) 

In the context of pair breaking, A should everywhere be 
replaced by the field-dependent spectroscopic gap il G . 

The gap dependence of Tr, t b , N th , and R 
can also be numerically evaluated directly using 
Eqs. dSTOl . ( ISTTb . dSBJ, and (|S9j without approxi- 
mations. The results are shown in FIG.IS7I confirming the 
linear relation between R and D, G , discussed in the main 
text and shown in FIG. 4c and 4d. 
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